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The thermal conductivity k of underdoped YBa2Cu30 H was measured in the T — > limit as 
a function of hole concentration p across the superconducting critical point at psc = 5.0 %. The 
evolution of bosonic and fermionic contributions to k was tracked as the doping level evolved contin- 
uously in each of our samples. For p ^ Psc, we observe a T 3 component in k which we attribute to 
the boson excitations of a phase with long-range spin or charge order. Fermionic transport, observed 
as a T- linear term in k which persists unaltered through psc, violates the Wiedemann- Franz law, 
since the electrical resistivity varies as log(l/T) and grows with decreasing p. 

PACS numbers: 74.25.Fy, 74.72. Bk, 72.15.Eb 



The nature of the ground state near the point where 
superconductivity vanishes in underdoped cuprates is of 
fundamental interest. While a pseudogap and strong 
phase fluctuations are well established properties of this 
state, a number of key issues remain open. Is there an 
order parameter? If so, what symmetry is broken? What 
are the elementary excitations? Theoretical models ad- 
dressing these questions involve, for example, the valence 
bond solid ci-density wave order stripe order 
or fluctuating vortex-antivortex excitations Recent 
measurements, on the other hand, show that magnetism 
can be present up to the point where T c vanishes J0, 0] , 
or that electrons can form a charge-ordered state [8| . 

In this Letter, we report a study of the evolution of 
heat and charge transport across the superconducting 
transition at psc in underdoped YB^CusOy. In this 
cuprate, the carrier concentration p in the CuC>2 planes 
is determined by the average oxygen content and the de- 
gree to which the dopant atoms are ordered in the chain 
layers |9[ . Oxygen ordering occurs at room temperature 
and the resulting continuous change of hole doping with 
time, called "time doping" offers a unique opportunity to 
track the evolution of one single specimen as a function 
of p without the uncertainty related to changes in cation 
disorder, or sample and contact characteristics (loj . 

The thermal conductivity k of any solid may be ex- 
pressed as the sum of contributions from fermions, bosons 
(other than phonons), and phonons, respectively: 



K b + K p = aT + f3T 3 + bT a 



(1) 



In the T — > limit where all inelastic scattering vanishes, 
the mean free path £ becomes independent of temper- 
ature and frequency and each contribution acquires its 
distinct asymptotic T dependence, reflecting that of the 
specific heat c(T) through k = ^cv£, where v is the car- 
rier velocity. One exception, as we shall discuss, is the 



phonon exponent a whose value typically falls below 3 
due to specular reflection at the sample boundaries. 

In practice, it can be difficult to distinguish Kb from 
k p unless one can be "tuned" while the other remains 
constant. This was shown in a recent study of un- 
doped Nd2Cu04 [ll[ where heat transport by ballistic 
magnons, revealed by a large T 3 contribution to k, could 
be "switched on" by field-tuning across a spin-flop transi- 
tion. Here, time doping is the tuning method and, like in 
Nd2Cu04, we have identified a sizable T 3 contribution to 
the thermal conductivity of YBa2Cri30 y that "switches 
on" when p is lowered below psc- We associate this term 
with non-phononic boson excitations characteristic of the 
non-superconducting phase. 

Heat and charge transport by fermions was also ob 



served. This follows the work of Sutherland et al. [12 1 
and Sun et al. 13]. In the former, a finite and doping- 
independent k / due to fermionic excitations was reported 
for p on both sides of psc- In contrast, for a y = 6.35 
sample Sun et al. found Kf = and a log(l/T) depen- 
dence of the electrical resistivity, suggesting an insulat- 
ing state. Here, we measured both properties on one 
single specimen as its doping evolves and found that K/ 
indeed remains finite and constant close to, and across, 
Psc- Since we also observe a log(l/T) electrical resis- 
tivity below psc, we conclude that at this doping the 
Wiedemann- Franz law is violated in YBa2Cu30 a . 

Ten high-purity single crystals of YBa 2 Cu30 y were 
measured, with oxygen content y = 6.0 (sample D), 
y = 6.31 (sample O), y — 6.32 (samples M, N), y = 6.33 
(samples J, K, L) and 6.35 (samples B, F, H) (for details, 
see [14j). Sample D is undoped and its hole concentration 
p is zero. In their oxygen-disordered state, the 6.33 sam- 
ples were not superconducting, but after 2 (29) days of 
time doping a T c of 0.1 (6.0) K was observed. For those 
and the 6.35 samples, p is determined using the known 



2 



dependence of T c on p 15| and may be read off Fig.[2j For 
example, p in sample L evolved in the sequence 4.7, 5.0, 
and 5.4 % after 0.25, 2, and 29 days, respectively. Since p 
varies logarithmically with time [9j , we can find the initial 
hole concentration of our 6.33 samples by extrapolating 
backward in time. As for the 6.31 and 6.32 samples, the 
lack of superconducting T c prohibits a precise determi- 
nation of p, whose initial value we loosely estimate to be 
near 4 %. T c is defined as the temperature where the 
resistivity becomes zero. The bulk T c is typically ~ 1 K 
lower, being the point where the susceptibility drop is 
complete. No sample was detwinned. 

Measurements were made using a four-contact AC 
method for electrical resistivity and a one-heater-two- 
thermometer technique for thermal conductivity [lrjj ]. 
Great care was taken to ensure a uniform current dis- 
tribution: very thin samples were used (4-10 /im thick, 
1000-2000 fim long) with contacts made by evaporating 
and diffusing gold pads on the sides of each sample. 

The thermal conductivity of samples L and B is shown 
in Fig. [U plotted as k/T vs T. One can identify three 
contributions: 1) a T 3 term, due to bosons, present be- 
low p = 5.3 % but not above, revealed by plotting the 
difference between successive anneals (upper inset); 2) a 
residual linear term kq/T, due to fermions, present in 
both samples; and 3) a phonon background which goes 
as T a with a < 3. We discuss each contribution in detail. 

Bosons. The data in Fig. [1] shows a generic feature of 
6.33 samples: K drops as p increases, by as much as 30 % 
at 0.6 K. When plotted vs T 2 (upper inset), the differ- 
ence k/T(4.7%) - k/T(5.4%) is linear below 450 mK and 
a fit to Ak/T = Aa + /3T 2 gives (3 = 476 ^W K~ 4 cm- 1 
and a negligible intercept (Aa = 0). A fit to the interme- 
diate difference, k/T (5.0%) - k/T(5.4%), yields a similar 
result, with (3 = 283 ^W K~ 4 cm" 1 and Aa = 0. Mea- 
surements on more underdoped samples, with y = 6.31 
and 6.32, show that this T 3 contribution to k persists at 
lower doping. In contrast, however, this term is no longer 
present above 5.2 %. This is shown in the lower inset of 
Fig. [TJ where sample B exhibits negligible change in k as 
its doping varies from 5.2 to 5.7 %. Collecting all sam- 
ples at all dopings in Fig. [2] reveals a systematic trend 
as p increases: (3 decreases roughly linearly until 5.2 %, 
at which point it stops evolving. We thus conclude that 
/3 = 0forp>5.2%. Note that the difference between 
5.0 %, where T c starts to grow, and 5.2 %, where [3 ex- 
trapolates to zero, is within the combined uncertainty in 
the bulk T c and the sampling density. 

A pure T 3 term in c(T) as T — > is the standard 
signature of bosonic excitations with finite dispersion in 
all three directions, i.e., as in any phase with long-range 
order. For this to be reflected in k(T), the mean free 
path needs to be independent of T, as in the ballistic 
regime where it is limited only by sample boundaries. 
In an anisotropic system like cuprates, one can expect 
k ~ LT 3 /v±v, where L is the sample size and v (v±) is 
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FIG. 1: Temperature dependence of k/T for three crystals 
of YBa2Cu30 y with different oxygen content y. Main panel: 
y — 6.33 (sample L) at p — 4.7 % (open circles) and 5.4 % 
(full circles), with a power-law fit to k/T — a + bT a ~ 1 (solid 
line), giving a — 2.0; y = 6.0 (p = 0) (full squares), with 
the same fit (dashed line), giving a — 2.4 (k/T is limited 
to T < 0.3 K and scaled by 1/4 for clarity). Upper inset: 
difference between k/T(4.7%) and k/T (5.4%), vs T 2 , with a 
linear fit. Lower inset: y = 6.35 (sample B) at p = 5.2 % 
(open triangles) and 5.7 % (full triangles), with a fit to the 
5.2% curve that yields a = 2.4. 



the in-plane (out-of-plane) acoustic boson velocity. Such 
ballistic acoustic bosons offer the most natural interpre- 
tation for the PT 3 term in k observed in YBa2Cu30j,. 

The appearance of a new boson mode below psc sug- 
gests that long range spin or charge order sets in where 
superconductivity vanishes. In the case of spin order, 
the /3T 3 term is ascribed to the magnons of an antifer- 
romagnetic state that develops below p$c ■ This scenario 
is motivated in part by the striking similarity with the 
well- understood antiferromagnet Nd2Cu04 [ll|: 1) the 
ballistic regime where Kb oc T 3 is nearly the same in 
both cuprates, up to around 0.5 K, and 2) the value of 
(3 is comparable when considering the difference in sam- 
ple size. Indeed, although (3 in Nd2Cu04 is 30 times 
larger than the value in YBa2Cu30 !/ at 4.7 %, the three 
6.33 samples are 10-15 times thinner. In further support 
of this scenario, recent neutron scattering measurements 
on YBa^CusOy show the in-plane spin correlation length 
growing rapidly as p — > p$c (from above) [13], and /itSR 
measurements show that the Neel temperature of low- 
doped YBa2Cu30 y declines rapidly as p — > 0.05 (from 
below) 0. In order to transport heat as T — > 0, we note 
that magnons would have to be gapless at the critical 
point. 

Another scenario is a transition to a Cooper pair 
Wigner crystal (PWC). Phase fluctuations of the d-wave 



3 




0.045 0.050 0.055 0.060 



Hole concentration 

FIG. 2: Boson and fermion contributions to the thermal 
conductivity of YBa2Cu30 H vs hole concentration p. Upper- 
panel: coefficient (3 of bosonic conductivity, Kb = [3T 3 . The 
full line and black dots are the resistive T c (p) (Note: the bulk 
T c lies somewhat lower). For p ^ 5.0 %, a finite resistive T c is 
observed, from which p is determined F° r V < 5.0 %, no 
T c is observed and we use a method described in the text to 
estimate p. Lower panel: coefficient a = kq /T of the fermionic 
conductivity, ft/ = aT. Error bars for f3 and kq/T reflect the 
statistical uncertainty in the fit. When comparing data from 
different samples, the uncertainty on the geometric factor, 
typically between 10 and 20 %, must be added. 



superconducting order parameter can lead to such order 
(lq ] and recent STM experiments on cuprates (e.g., Q) 
have been interpreted as supporting evidence. Pereg- 
Barnea and Franz [lj| recently showed that a PWC has 
gapless transverse modes of vibration that contribute a 
T 3 term to k, with a coefficient they estimate to be in 
broad agreement with the experimental value of f3 re- 
ported here. We note that the fermion sector of the same 
theory ([3, [T^ and references therein) is qualitatively 
consistent with other aspects of our study, as we shall 
see below, in the sense that charge is localized in the 
PWC, but spin can remain mobile, in the form of neutral 
fermionic spin excitations with a Dirac spectrum. 

Fermions. The fermionic contribution to k is obtained 
from the T = intercept kq/T. For sample L at 5.4 %, 



a power-law fit to k/T = a + bT a 1 up to 0.6 K yields 
a = kq/T = 40 ± 1 /xW K~ 2 cm" 1 , with a = 2.0. As 
already mentioned, the difference An/T between k/T at 
4.7 % and k/T at either 5.0 or 5.4 % yields a negligible 
intercept, showing that kq/T for L does not change with 
doping upon crossing psc- I n Fig- El this and the data 
for five other samples are shown, obtained using the same 
procedure: a power-law fit for p ^ 5.3 % and a T 2 fit to 
the difference for p < 5.3 %. It is remarkable that all 
six samples measured at p ~ 5.3% gave values within 
10 % of each other: k /T = 38 ± 3 pW K~ 2 cur 1 . 
This highly reproducible finding confirms the presence 
of fermionic excitations below psc, as reported recently 
fl3 |. A magnetic field up to 10 T has negligible effect on 
k for all samples and dopings. 

The fact that kq/T is nearly constant across psc sug- 
gests that the spectrum of nodal quasiparticles associated 
with the (i-wave superconducting gap does not change 
when superconducting order vanishes. This is what one 
might expect of a "nodal metal" , a state characterized in 
recent ang le-resolved photoemission spectroscopy experi- 
ments [20| . Measurements on samples with y = 6.31 and 
6.32 show that k /T eventually becomes negligible as p 
approaches 4 %, so that fairly soon this nodal metal phase 
either comes to an end (e.g., through the appearance of 
a gap at the nodes) or its fermions become localized. 

Our observation of a crossover in kq /T, from vanishing 
near 4 % to non- vanishing at 4.7 % (or below), is quali- 
tatively consistent with the results of Sun et al. [1 31 ] , who 
put the crossover between their more underdoped sam- 
ple, with p ~ 4.5 %, and the next one up, with p ~ 5.5 % 
(T c = 7 K). They conclude that the crossover occurs 
halfway, i.e., at 5 %. With the greater resolution con- 
ferred by time doping a single sample across Psc, w e 
can show unambiguously, as already done in [12] . that 
there is a region of the phase diagram with delocalized 
fermions and no superconductivity. Note that beyond 
this issue, there appears to be some discrepancy between 
their data and ours. While both yield the same value of 
kq/T at p ~ 6 %, theirs gives a significantly lower value 
at p ~ 5.5 %. We can only speculate that this "metal- 
insulator" crossover is pushed to slightly higher p in their 
samples because of a higher impurity level (coming from 
growth in zirconia vs barium-zirconate crucibles) . 

Our observation of fermionic heat conduction begs the 
question as to what happens to charge conduction. Fig. [3] 
shows the normal state electrical resistivity p of sample 
L at 4.7 and 5.0 % in H = and 10 T, respectively. We 
see that the electrical resistivity drops with increasing 
doping (at all T), as expected from the increase in carrier 
density, and that it exhibits a perfect log(l/T) behaviour 
(see inset), over two decades in T. The latter shows a 
clean (and strong) yet unusual tendency towards charge 
localization, as already reported by Sun et al. [13J . 

The measured electrical resistivity is in sharp contrast 
with the corresponding fermionic heat transport, plot- 
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FIG. 3: Electrical resistivity p of non-superconducting 
YBa2Cu3 06.33 at p = 4.7 (open triangles) and 5.0 % (full 
triangles) in H = and 10 T, respectively. The correspond- 
ing fermionic thermal resistivity is plotted in the same units, 
as LqT/kq (open and closed circles). Inset: p(T) at p — 4.7 % 
(H — T) on a logarithmic scale, with linear fit. 



ted in Fig. [3] in units of resistivity, i.e., p — LqT/kq, 
where Lq = (7r 2 /3)(kB/e) 2 . The significant difference 
in magnitude reveals an unambiguous violation of the 
Wiedemann-Franz law (WFL), k/T = L /p as T -> 0. 
Beyond the quantitative discrepancy, it is the qualitative 
trend that is compelling: in the same sample, p changes 
with p while kq/T remains constant. This is consistent 
with nodal (and neutral) fermions governed by a Dirac 
spectrum unrelated to the carrier density, much like d- 
wave quasiparticles in the superconducting state (l6| . 

In the other two instances of WFL violation, ob- 
served in Pr2_ s Ce cc Cu04_,5 (PCCO) at optimal doping 
(x = 0.15) [2l[ (see (22[) and in slightly underdoped 
Bi 2+a: Sr2- :c Cu06+5 (Bi-2201) H, k /T was also found 
to exceed Lq/ p. Given the widely different magnitudes 
and T-dependencies of the electrical resistivity in the 
three cases, more work is needed to elucidate the cause 
of the violation, and determine whether it is the prox- 
imity to a metal-insulator crossover that is relevant, as 
proposed by Proust et al. for Bi-2201 [23], or rather the 
proximity to a magnetic quantum critical point, such as 
that reported for PCCO at x = 0.16 [II]. 

Phonons. Unlike the fermion and boson parts, k/ and 
Kb, whose magnitudes are highly reproducible for the six 
samples, the phonon background k p — n — Kf — k& = bT a 
is sample dependent (but constant as a function of p for 
each specimen), with a ranging from 1.8 to 2.5. This 
is consistent with k p being dominated by boundary scat- 
tering and specular reflection, both dependent on surface 
dimensions and quality [16j | . Indeed, k p in cuprates never 



bT a , with 2 < a < 3, 
The most convincing illustration of this sub-T 3 depen- 
dence comes from data on undoped cuprates, insulators 
in which phonons are the only carriers of heat (since 
magnons are gapped). In Fig. [1] we show the thermal 
conductivity data for a typical crystal of YBa2Cu306.o 
(undoped YBa2Cu30 y ) along with the corresponding fit 
to k/T = a + bT a ~ x which describes the data well from 
50 to 500 mK, and yields a ~ and a = 2.4. 

In summary, by carefully tuning the hole concen- 
tration p of underdoped YBa2Cu30 y crystals in- situ 
through the superconducting critical point at psc, we 
have identified the low-lying excitations that characterize 
the superconductor-insulator transition, from heat trans- 
port at T — > 0. A novel boson excitation is seen to emerge 
below psci most likely the boson of a phase with long- 
range spin or charge order. Possible interpretations in- 
clude the magnons of an antiferromagnetic phase or the 
vibrational modes of a Cooper pair Wigner crystal. On 
the other hand, the fermionic component of heat conduc- 
tion persists unchanged through psCi even as the charge 
conduction decreases and localizes. This suggests that 
fermion excitations with nodal (and neutral) character 
survive outside the superconducting phase, at least in the 
vicinity of the critical point. It is interesting to consider 
the evolution of charge transport below psc in relation 
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to the superfluid density above psc 
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